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Abstract 



By considering the long-wave limit of the regularized long wave (RLW) equa- 



tion, we study its multiple-time higher-order evolution equations. As a first 
result, the equations of the Korteweg-de Vries hierarchy are shown to play a 



^ \ crucial role in providing a secularity-free perturbation theory in the specific 

case of a solitary-wave solution. Then, as a consequence, we show that the 

related perturbative series can be summed and gives exactly the solitary-wave 

solution of the RLW equation. Finally, some comments and considerations are 

made on the N-soliton solution, as well as on the limitations of applicability 

of the multiple scale method in obtaining uniform perturbative series. 
03.40.Kf 
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I. INTRODUCTION 



The regularized long wave (RLW) equation, 



u t + u x - u xxt - 6uu x = , 



(1) 



also known as Peregrine ]I[ or Benjamin-Bona- Mahony equation, was originally intro- 
duced as an alternative for the Korteweg-de Vries (KdV) equation 



u t + u x + u. 



xxx 



— 6uu 



x 



0. 



(2) 



Despite having quite different dispersion properties, these two equations present some inti- 
mate relationship. For example, the linear dispersion relation of the RLW equation is 



which, by the way, is the same as the dispersion relation of the shallow water wave equa- 
tion ||. For long waves, k is small and u(k) can be expanded according to 



The first two terms of this expansion coincide exactly with the complete linear dispersion 
relation of the KdV equation. Thus, for sufficiently long waves, the traveling-wave solutions 
of Eqs. (HD and @ are expected to be quite similar. Despite of this, there is a deep difference 
between these two cases since a polynomial is definitely not equivalent to an infinite series. 
This difference, which appears when higher order terms of the dispersion relation expansion 
are considered, might also show up at higher order approximations in a perturbation theory. 

As is well known, the KdV equation appears as governing the first relevant order of an 
asymptotic perturbation expansion describing weakly nonlinear dispersive waves. However, 
to make sense of it as really governing such waves, the large time behaviour of the perturba- 
tive series must be analyzed [f| . In other words, one needs to study the evolution equations 
of the higher order terms of the perturbative expansion to check for the existence or not of 



uj{k) 



k 



(3) 



1 + k 2 




(4) 
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secular producing terms. This study, usually neglected in the derivation of the KdV equa- 
tion, is essential to guarantee the uniformity of the perturbative expansion, rendering thus 
a real meaning to the KdV equation. 

Motivated by the above considerations, we are going in this paper to apply a multiple 
time version []5J of the reductive perturbation method to study long waves as governed by 
the RLW equation. As we are going to see, the KdV equation appears at the lowest relevant 
order of the perturbative scheme. Then, by assuming a solitary-wave solution for the KdV 
equation, we consider the higher order approximations and we show that the corresponding 
solitary-wave related secular producing terms can be eliminated from every order of the 
perturbative scheme. The equations of the KdV hierarchy, which appear as a consequence of 
natural compatibility conditions, are shown to play a crucial role in the process of eliminating 
the secular producing terms. Once a secularity-free perturbative series is obtained, we show 
that it may be summed to give the exact solitary-wave solution of the RLW equation. We 
then close the paper with a discussion on the N-soliton well as on the limitations of 

the multiple scale method. 

II. MULTIPLE-TIME FORMALISM 

To study the long-wave limit of the RLW equation we put 

k = en , (5) 

with e a small parameter. In this limit, the dispersion relation (§) can be expanded as 

uj(k) = en- e 3 K 3 + e 5 n 5 - e 7 n 7 H . (6) 

Accordingly, the solution of the corresponding linear RLW equation can be written in the 
form 

u = a exp i[kx — u(k)t] = a exp i kc(x —t) + e 3 K 3 t — e 5 K 5 t + e 7 K 7 t + ■ ■ • , (7) 
where a is a constant. As given by this solution, we define now a slow space 



£ = e(x - t) , (8) 
as well as an infinity of properly normalized slow time variables: 

r 3 = e 3 t ; r 5 = -e 5 t ; r 7 = e 7 t ; etc. (9) 



Consequently, we have 



1 = 4- (10) 



and 



— = -e he 3 e he 7 . (11) 

It is important to notice that the introduction of slow time variables normalized according 
to the dispersion relation expansion are such that they allow for an automatic elimination 
of the solitary-wave related secular-producing terms appearing in the evolution equations 
for the higher order terms of the wave-field ||. 

III. PERTURBATION THEORY 

The perturbative scheme consists of making the expansion 

u = e 2 u = e 2 (u + e 2 u 2 + e 4 u 4 + •••), (12) 

and substituting it, together with Eqs.([T0|) and (|TTD, into the RLW equation ([[]). The result 
is the multiple-time equation 



<9t 3 9t 5 J <9£ 2 \ <9£ <9r 3 ' dri 
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We proceed then to an order-by-order analysis of this equation. 
At the lowest order, we obtain 



3e 3 ^-{u 2 + 2e 2 u u 2 + ---) =0. (13) 



«0r 3 = F 3 = -«0& + 6U U £ = , (14) 
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which is the KdV equation. Introducing an operator L, whose action on any component u n 
is given by the linearized KdV operator 

Lu n = u nT3 + Un^ - G(u U n )(: , (15) 

the KdV equation ( |T4"D can be rewritten in the form 

Luq = —6uoUo£ . (16) 

Our interest in this paper is concerned to solitary-waves. Thus we assume uq to be the 
solitary- wave solution of the KdV equation (|14D , 

uo = -2 k 2 sech 2 # 3 , (17) 

where #3 = «[£ — 4k 2 t 3 ]. In this case, Eq. (|16|) becomes 

Lu = 48 k 5 sech 4 6 3 tanh 6 3 . (18) 

In the next order, we obtain the equation 

Lu 2 = u 0rM + u 0T5 . (19) 

The evolution of uq in the time T3 is given by the KdV equation (|i"4j) , but the evolution 
of Uo in the time t$ is not known up to this point. However, the multiple-time formalism 
introduces constraints which determine uniquely the evolution of in any higher-order 
time 0. To see how this is possible, let us make the following considerations. 

First, to have a well ordered perturbative scheme, we impose that each one of the equa- 
tions describing the higher-order times evolution of uq be e-independent when passing from 
the slow («, Mo, £, T 2n+ i) to the laboratory coordinates (k, u, x, t). This will select all possible 
terms to appear in UQ T2n+1 . For instance, the evolution of Uq in the time T5 is restricted to 
be of the form 

uot 5 = "5^0(50 + fcuouo^ + 75M05M05S + 5 5 ulu z , (20) 
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where 0*5, ^5, 75 and £5 are constants to be determined. Then, by imposing the natural (in 
the multiple time formalism) compatibility condition 

(«0r 3 ) T5 = (^0r 5 ) r3 , (21) 

it is possible to determine the above constants in terms of a§, which is left as a free- 
parameter. As it can be verified ||, the resulting equation is the 5th order equation of the 
KdV hierarchy: 

u 0t5 = F 5 = M 0(5 £) - 10m Mo£& - 20m £Mo& + 30u 2 u 0i . (22) 

The right-hand side of this equation would in principle appear multiplied by the free param- 
eter a 5 , which would account for different possible normalizations of the time r 5 . However, 
since we have already defined the slow time normalizations, this parameter was taken to be 
1 in order to have an agreement with the normalizations introduced in Eq.([|). This is an 
important point since, as we have already said, it allows for an automatic elimination of the 
solitary- wave related secular producing terms appearing in the right-hand side of Eq.(|T^). 
These terms, when uq is assumed to be a solitary-wave of the KdV equation, are always of 
the form J7J 

U [(2n+m; 71 = 0,1,2,... . (23) 

Thus, using Eqs. (|T4"D and (^) to describe respectively u 0t3 and m 0t5 , Eq.(|T9|) becomes 

Lu 2 = -2u £U tt ~ 4wo«o&£ + 30w 2 «o£ • (24) 

We notice in passing that the substitution of Eqs. (|i~4[) and (0), respectively for m 0t3 and 
Uq T5 , with the properly normalized slow times allowed for an automatic elimination of all 
solitary- wave related secular producing terms of Eq.([L^). In fact, Eq.(|2~4|) does not present 
any secular-producing term anymore. Moreover, we see that at this order uq must satisfy 
simultaneously the first two equations of the KdV hierarchy, respectively in the slow-times 
T3 and T5. Introducing the general definition 
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#2n+i = « [e - 4^ 2 r 3 + 16/£ 4 r 5 - • • ■ + (-l) n (2«) 2n r 2n+1 ] , (25) 
such a solitary- wave is given by 

m = -2 sech 2 # 5 , (26) 

and Eq.(^) becomes 

Lu 2 = 192 k 7 sech 4 6 5 tanh 6> 5 . (27) 

Assuming a vanishing solution for the associated homogeneous equation, we can write the 
solution of this equation in the form 

u 2 = 4k 2 m , (28) 

with uq given by fl26|) . 

We proceed then to the next order, where we get 

Lu A = -u 0t7 - u 0T5 ^ + u 2T5 + u 2t3 ^ + 6u 2 u 2 {: . (29) 

Following the same scheme used before, we can use the compatibility condition 

(u 0t3 ) t7 = (u 0t7 ) t3 (30) 

to obtain the evolution of uq in the time Tj. It is given by 

u 0r7 = F 7 = - u 0(7O + 14w Wo(5£) + 42n ^ 0(4C) + U0(v ) 3 Vqz 

+ 70u ^u ^ - 280u o u O £U tt ~ 70 ( w o?) 3 - 70u 2 u ^ , (31) 

which is exactly the 7th order equation of the KdV hierarchy. At this order, therefore, the 
solitary-wave must satisfy simultaneously the first three equations of the KdV hierarchy, 
respectively in the times T3, T5 and r 7 . This means that 

u = -2K 2 sech 2 9 7 . (32) 
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Now, by using Eq. (|28|) to express u 2 , and the equations of the KdV hierarchy to express u 0tt , 
u 0t5 and u 0t3 , all secular producing terms of Eq.(^9|) are automatically eliminated. Then, 
substituting the solution (|32|), Eq.([29|) becomes 

Lu 4 = 768 k 9 sech 4 6 7 tanh 6 7 . (33) 

Again, by assuming a vanishing solution for the associated homogeneous equation, the so- 
lution of this equation can be written as 

Ui = (Ak 2 ) 2 u , (34) 

with Uq given now by Eq. fl32|) . 

This procedure can be repeated up to any higher order. In other words, we can use the 
compatibility condition 

(%r 3 ) = (u r 2n+1 ) (35) 

to obtain the evolution of u in the time T2 n +i, which will turn out to be the (2n + l)th 
equation of the KdV hierarchy. In this case, Uq will represent a solitary-wave satisfying 
simultaneously the first n equations of the KdV hierarchy: 

Uq = —2k 2 Sech 2 6*2n+l • (36) 

The resulting secularity-free evolution equation at this order will be 

L u 2n = 3 (4)™+ 2 (n) 2n+5 sech 4 9 2n+3 tanh 9 2n+3 . (37) 

Assuming a vanishing solution for the associated homogeneous equation, the solution to this 
equation can be written in the form 

u 2n = (4K 2 ) n u , (38) 

with uq given by Eq. (p6|) . Extending this procedure ad infinitum, u will represent a solitary- 
wave satisfying simultaneously all equations of the KdV hierarchy, and we obtain an exact 
solution for the RLW equation. 
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IV. RETURNING TO THE LABORATORY COORDINATES 



Let us take the solutions u 2n and substitute them in the expansion flT2|). Putting u in 
evidence, we get 



u = e u 



1 + 4eV + 16e 4 « 4 + 6Ae 6 K 6 + 



(39) 



Now, the above series can be summed: 

1 + 4eV + 16eV + 64e 6 K 6 + 

Therefore, we get the RLW exact solution 

2e 2 K 2 



1 - Ae 2 K 2 



(40) 



u 



sech 



- 4/t d r 3 + 16k°t 5 - 64k't 7 + 



(41) 



1 - Ae 2 K 2 

Then, by using Eqs.@, (§) and (|), we can rewrite u in terms of the laboratory coordinates 
(k,x,t). The result is 

2k 2 



u 



l-Ak 2 



sech z 



kx — k(l + Ak + 16F + 64A; b + • • • )t 



(42) 



Using again Eq.(^0|) with en = k, we get finally 

t 



u = —a sech 



k x 



l-Ak 2 



2k 2 



l-Ak 2 



(43) 



which is the solitary- wave solution of the RLW equation (|l]). 
The RLW equation has another solution, given by 



u = b tanh 



k ( x 



l + 8k 2 



2k 2 



l + 8k 2 



(44) 



In fact, it is easy to see that Eq.(P is invariant under the transformation 

t' = a~ t ; x' = x ; u' = b — au , (45) 

where, if u is given by Eq. fl43|) , u' turns out to be the solution given by Eq.([H]). By following 
the same procedure used to obtain the RLW solitary- wave solution ([43]), it is also possible 
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to use the multiple-time perturbative scheme to obtain the solution (fHj). This is done by 
choosing 



u 



2/t 2 tanh 2 (k£ - 8k 1 



T 3 



(46) 



instead of ( |TTD as the solution for the KdV equation (13). As higher orders are reached, 
this uq is required to satisfy also the higher order equations of the KdV hierarchy, which 
amounts to include dependences on the higher order times r 5 , r 7 , etc. However, there is an 
important difference: the secular-producing term in each order of the perturbative scheme 
will come not only from the linear term, but from both the linear and the nonlinear terms. 
As a consequence, the slow time normalizations obtained from the linear dispersion relation 
expansion will not be able to remove the secular-producing terms in this case. In other words, 
new slow time normalizations will be needed to get a secularity-free perturbative series. 
These new normalizations can be easily found by properly choosing the free-parameters left 
at each order of the perturbation scheme ||. After doing that, we obtain the following 
perturbative series for u: 



u 



8k 2 + 64k 4 



tanh 



kx - k (l - 8A; 2 + 64k 4 



(47) 



Like in the previous case, these series can be summed, resulting in 



2k 2 



u 



Ik 2 



tanh 



k I x 



t 



8k 2 



(48) 



which is the solution ( fPj) of the RLW equation. As already said, however, a new slow time 
normalization is needed in this case to get a secularity-free perturbative series, which is 
different from that obtained from the dispersion relation expansion. 



V. STUDY ON THE APPLICABILITY OF THE MULTIPLE SCALE METHOD 

The multiple scale method is not always able to remove all the secular-producing terms of 
a perturbative series |5J. In some cases, nonintegrable effects may preclude the existence of 
uniform asymptotic expansions. Considering that the RLW is nonintegrable, the purpose of 
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this section will be to make a brief discussion on how those effects appear in the higher order 
terms of the perturbative series for the specific case of the RLW equation. The approach we 
are going to use is that developed by Kodama and Mikhailov ||. 
Let us start by defining slow variables according to 



u = ev ; £ = e L/2 (x-t) ; r 3 = e i/2 t. 
In these new coordinates, and up to terms of order e 2 , Eq.(0) becomes 



v T3 = 3v 2 — V££ + ed& (3v 2 — v^A + e 3 <9(4^) (3v* 
Then, we make a near identity transformation || given by 



v = w + e$(u>) + e 2 ^(w) + 0(e 



(49) 



(50) 



(51) 



where, by reasons of scaling- weight invariance, the differential polinomials $ and ty, which 
are allowed to be nonlocal, can involve only the following terms: 



$ = aw 2 + /3w££ + ■yw^d 1 w , 



(52) 



\& = aw 3 + b(w^) 2 + cww££ + dw(4£) + eww^d w 

+ fw^d~ 1 (w 2 ) + gw^d~ x w + hw^d^w) 2 



(53) 



Substituting into (50), we obtain 



w T3 = K 3 + eK 5 + e 2 K 7 + 



(54) 



with 



K 3 = <%M , 



K 5 = (M 1 + %M ) - ^ (d c M ) 



K 7 = dAM 2 + %Mi + d m M ) - 



5^ 

5w 

Sw 



(%M ) 

a 5 {m x + %m ) 



5$ 



(55) 
(56) 



(57) 
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where we have introduced the notation: 

M = 3w 2 - , (58) 

M l = 6w$ - % , (59) 

M 2 = 3$ 2 + 6wV - % . (60) 

At order e° we find 

K 3 = F 3 = 6ww^ - , (61) 

that is, K 3 is the symmetry of order e° of the KdV equation. At the next order, by properly 
choosing a, j3 and 7, we find 

K 5 = F 5 , (62) 

with F 5 defined by Eq.(kdv5).This means that there exists a near-identity transformation 
(pT|)-([5B|) such that K5 is the symmetry of order e of the KdV equation. In the first two 
orders, therefore, no problems appear. This is a general result that holds for any equation, 
not only for the particular case of the RLW equation. It is in the next order that the so 
called obstacles || show up. In fact, in the next order we get 

K 7 = F 7 + 0{w) , (63) 

with F 7 defined by Eq. (|31|) , and 0{w) representing the obstacle, which is given by 

°{ w ) = ( -"V ~ 3 #) ww m + (--§■ - 3c - 24d - 3g) w^w {40 





( 32 


= I 
+ i 


v~y ~ 

( 508 




I 3 " 


+ i 


f700 
V 3 


+ i 


/158 

(3 1 



3 

f 6a + 2f - 18c/) w 3 W£ + (22 - 3c - 6b - 60d) w^w^ 

\ /224 \ 

18a - 12c - 6f + 72g) ww^w^ + I — ~ 3/ + 21g\ w 2 w^ 



6a - 6b - 3f + 18g ) {wrf . (64) 
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The important point is that, for an arbitrary KdV hierarchy solution w, it is not possible 
to choose a,b, . . . ,g in such a way to have a vanishing obstacle. However, as an explicit 
calculation easily shows, when w is a solitary-wave solution of the KdV hierarchy, there is a 
near-identity transformation leading to 0(w) = 0. 

The above considerations are important in the sense that they clarify the results obtained 
in the previous sections concerning the solitary- wave related secularities. But, at the same 
time, they put in evidence the limitations of the perturbative scheme which, as we now 
know, can not be extended to the two-or-more soliton solutions in the non-integrable case. 
On the other hand, for integrable systems, like for example the shallow water wave equation, 
the multiple scale method will be able to handle both, the solitary-wave and the N-soliton 



related secularities [11] since no obstacles will be present in either case. 



VI. FINAL REMARKS 

We have applied a multiple-time version of the reductive perturbation method to study 
the solitary- wave solution of the RLW equation. As it has already been shown || , the use of 
multiple time-scales allows for the elimination of all solitary-wave related secular-producing 
terms appearing in the evolution equations of the higher-order terms of the wave-field. More- 
over, it has also been shown || that these secularities are automatically removed if the slow 
time-scales are normalized according to the long-wave expansion of the dispersion relation of 
the original equation. By using this strategy, we have succeeded in expressing the solitary- 
wave solution of the RLW equation as a sum of solitary-waves satisfying simultaneously, in 
the slow coordinates, all equations of the KdV hierarchy. Similar results have been shown 



to hold also for the Boussinesq Jl0| and the shallow water wave |TTJ equations. However, 
while in these two cases the solitary-wave solution was obtained due to a truncation of 
the perturbative series, the RLW solitary-wave was obtained by summing the perturbative 
series. 

To finish, let us make the following considerations. If we assume the RLW equation 
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to be an exact model equation, as we have in fact done, the KdV equation appears as its 
long wave leading order approximation. This is one more confirmation of the widely known 
property of the KdV equation, which states that it holds a unique, privileged and universal 
meaning in the sense it appears as the leading order approximation of any weakly nonlinear 
dispersive systems, as for example that represented by the RLW equation. From this point 
of view, the old dispute |12| on the equivalence of the RLW and the KdV equations would 
be made on a different ground since the RLW equation should be compared not to the 
KdV equation, but to the whole set of equations of the KdV hierarchy. In other words, the 
RLW equation should be compared not to its leading order approximation, but to the whole 
perturbative series. And according to our results, as far as solitary-waves are concerned, 
the RLW equation is indeed equivalent to the KdV hierarchy since a solitary-wave of the 
RLW equation is nothing but an infinite series given by the sum of solitary-waves satisfying 
simultaneously all equations of the KdV hierarchy, each one in a different slow time variable. 
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